石英晶体板非线性高频振动的拓展伽辽金法分析

针对考虑几何和材料非线性的石英晶体板厚度剪切振动和弯曲振动的方程组,利用扩展伽辽金法对该方程组进行转化和求解,分别获得了强烈耦合的厚度剪切振动模态和弯曲振动模态的频率响应关系,绘制了不同振幅比和不同驱动电压影响下的频率响应曲线图。数值计算结果表明可以选取石英晶片的最佳长厚比尺寸来避免两种模态的强烈耦合。驱动电压的变化将引起石英晶体谐振器厚度剪切振动频率的明显改变,必须将振动频率的漂移值控制在常用压电声波器件的允许值之内。扩展伽辽金法对石英晶体板非线性振动方程组的求解为非线性有限元分析和偏场效应分析奠定了基础。     

Analysis of a new multispecies tumor growth model coupling 3D phase-fields with a 1D vascular network

In this work, we present and analyze a mathematical model for tumor growth incorporating ECM erosion, interstitial flow, and the effect of vascular flow and nutrient transport. The model is of phase-field or diffused-interface type in which multiple phases of cell species and other constituents are separated by smooth evolving interfaces. The model involves a mesoscale version of Darcy’s law to capture the flow mechanism in the tissue matrix. Modeling flow and transport processes in the vasculature supplying the healthy and cancerous tissue, one-dimensional (1D) equations are considered. Since the models governing the transport and flow processes are defined together with cell species models on a three-dimensional (3D) domain, we obtain a 3D–1D coupled model.     

Subset representations and eigenvalues of the universal intertwining matrix

We solve a combinatorial question concerning eigenvalues of the universal intertwining endomorphism of a subset representation.     

A simple three‐wave approximate Riemann solver for the Saint‐Venant‐Exner equations

Erosion and sediments transport processes have a great impact on industrial structures and on water quality. Despite its limitations, the Saint‐Venant‐Exner system is still (and for sure for some years) widely used in industrial codes to model the bedload sediment transport. In practice, its numerical resolution is mostly handled by a splitting technique that allows a weak coupling between hydraulic and morphodynamic distinct softwares but may suffer from important stability issues. In recent works, many authors proposed alternative methods based on a strong coupling that cure this problem but are not so trivial to implement in an industrial context. In this work, we then pursue 2 objectives. First, we propose a very simple scheme based on an approximate Riemann solver, respecting the strong coupling framework, and we demonstrate its stability and accuracy through a number of numerical test cases. However, second, we reinterpret our scheme as a splitting technique and we extend the purpose to propose what should be the minimal coupling that ensures the stability of the global numerical process in industrial codes, at least, when dealing with collocated finite volume method. The resulting splitting method is, up to our knowledge, the only one for which stability properties are fully demonstrated.     

Local discontinuous Galerkin method for a nonlocal viscous conservation laws

The purpose of this article is to investigate high‐order numerical approximations of scalar conservation laws with nonlocal viscous term. The viscous term is given in the form of convolution in space variable. With the help of the characteristic of viscous term, we design a semidiscrete local discontinuous Galerkin (LDG) method to solve the nonlocal model. We prove stability and convergence of semidiscrete LDG method in L² norm. The theoretical analysis reveals that the present numerical scheme is stable with optimal convergence order for the linear case, and it is stable with sub‐optimal convergence order for nonlinear case. To demonstrate the validity and accuracy of our scheme, we test the Burgers equation with two typical nonlocal fractional viscous terms. The numerical results show the convergence order accuracy in space for both linear and nonlinear cases. Some numerical simulations are provided to show the robustness and effectiveness of the present numerical scheme.     

Generalized Jacobi spectral Galerkin method for fractional pantograph differential equation

This work is concerned with the extension of the Jacobi spectral Galerkin method to a class of nonlinear fractional pantograph differential equations. First, the fractional differential equation is converted to a nonlinear Volterra integral equation with weakly singular kernel. Second, we analyze the existence and uniqueness of solutions for the obtained integral equation. Then, the Galerkin method is used for solving the equivalent integral equation. The error estimates for the proposed method are also investigated. Finally, illustrative examples are presented to confirm our theoretical analysis.     

Novel local smoothness indicators for improving the third‐order WENO scheme

The local smoothness indicators play an important role in the performance of a weighted essentially nonoscillatory (WENO) scheme. Due to having only 2 points available on each substencil, the local smoothness indicators calculated by conventional methods make the third‐order WENO scheme too dissipative. In this paper, we propose a different method to calculate the indicators by using all the 3 points on the global stencil of the third‐order WENO scheme. The numerical results demonstrate that the WENO scheme with the new indicators has less dissipation and better resolution than the conventional third‐order WENO scheme of Jiang and Shu for both smooth and discontinuous solutions.     

Domain decomposition methods in a geometrical multiscale domain using finite volume schemes

The heat equation is solved by using a finite volume discretization in a domain that consists of a two-dimensional central node and several one-dimensional outgoing branches. Several interface connection options to match the submodels set on the node and on the branches, with or without continuity, are looked at. For each of them, a monolithic scheme is defined, and existence and uniqueness of the solution is proved. New schemes are deduced, which are obtained through domain decomposition methods in the form of interface systems, with one or two unknowns per interface. A comparative systematic study is carried out from an algebraic and numerical point of view according to the interface conditions: Dirichlet, Neumann, or Robin. An efficient diagonal preconditioning is proposed.     

Accurate iteration-free mixed-stabilised formulation for laminar incompressible Navier–Stokes: Applications to fluid–structure interaction

Stabilised mixed velocity–pressure formulations are one of the widely-used finite element schemes for computing the numerical solutions of laminar incompressible Navier–Stokes. In these formulations, the Newton–Raphson scheme is employed to solve the nonlinearity in the convection term. One fundamental issue with this approach is the computational cost incurred in the Newton–Raphson iterations at every load/time step. In this paper, we present an iteration-free mixed finite element formulation for incompressible Navier–Stokes that preserves second-order temporal accuracy of the generalised-alpha and related schemes for both velocity and pressure fields. First, we demonstrate the second-order temporal accuracy using numerical convergence studies for an example with a manufactured solution. Later, we assess the accuracy and the computational benefits of the proposed scheme by studying the benchmark example of flow past a fixed circular cylinder. Towards showcasing the applicability of the proposed technique in a wider context, the inf–sup stable P2–P1 pair for the formulation without stabilisation is also considered. Finally, the resulting benefits of using the proposed scheme for fluid–structure interaction problems are illustrated using two benchmark examples in fluid-flexible structure interaction.